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We investigate the effect of anisotropy and weak dipolar interactions on the mag- 
netization of an assembly of nanoparticles with distributed magnetic moments, i.e., 
assembly of magnetic nanoparticles in the one-spin approximation, with textured 
or random anisotropy. The magnetization of a free particle is obtained either by a 
numerical calculation of the partition function or analytically in the low and high 
field regimes, using perturbation theory and the steepest-descent approximation, re- 
spectively. The magnetization of an interacting assembly is computed analytically 
in the range of low and high field, and numerically using the Monte Carlo technique. 

Approximate analytical expressions for the assembly magnetization are provided 
which take account of the dipolar interactions, temperature, magnetic field, and 
anisotropy. The effect of anisotropy and dipolar interactions are discussed and the 
deviations from the Langevin law they entail are investigated, and illustrated for 
realistic assemblies with the lognormal moment distribution. 

PACS numbers: 75.50.Tt - 75.10.Hk - 05.20.y 

I. INTRODUCTION 

In order to investigate the fundamental properties of magnetic nanoparticles and the 
novelties they exhibit, new materials had to be made and characterized [see, e.g., the recent 
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review articles . Nowadays, there are mainly two prototypes of nanoparticle samples: i) 
assemblies of, e.g., cobalt, Nickel or maghemite nanoparticles |3| embedded in a non-magnetic 
matrix with volume distribution and randomly oriented easy axes, with negligible-to-strong 
dipole-dipole inter-particle interactions (DDI), depending on concentration; ii) isolated sin- 
gle particles of cobalt or nickel measured by the technique of /z-SQUID 4]. Technological 
applications require to some extent ever denser assemblies and thus smaller particles. How- 
ever, this leads to a dilemma because small particles become superparamagnetic at even low 
temperatures, and an optimum material [with appropriate anisotropy and other physical 
parameters] has still to be devised. Moreover, high density entails strong DDI among the 
particles, and in technological applications such as magnetic recording, this is an issue of 
special importance because DDI have been widely recognized as being responsible for the 
deterioration of the signal-to-noise ratio [see e.g., Refs. 0,0 and references therein]. 

Experimentally, investigation of the effect of dipolar interactions in nanoparticle assem- 
blies has revealed many new phenomena pertaining to the collective behavior of the particles, 



notably the so-called spin-glass state at low temperature in concentrated assemblies 
owing to the long-range of inter-particle DDI. It has also been observed that the field behav- 
ior of the temperature T max at the maximum of the zero-field-cooled magnetization strongly 
depends on the concentration of the assembly More precisely, the maximum 

of ^max as a function of the applied field observed in dilute samples disappears when the 
concentration of the latter is increased. Today, there arises a more fundamental issue about 
assemblies of nanoparticles that concerns the understanding of the interplay between the in- 
trinsic properties, such as those pertaining to surface effects, and extrinsic or collective effects 
stemming from the long-range DDI. Many research groups have experimentally studied this 
interplay in cobalt and maghemite particle assemblies. Measurements of the magnetization 
at high fields performed on the 7-Fe2C>3 nanoparticles 

[see Q for cobalt particles] 

have shown that the magnetization is strongly influenced by surface effects, depending on 
the particle size. For instance, Fig. 1 of Ref. shows that i) there is a sudden increase 
of the magnetization as a function of the applied field when the temperature reaches 70 K, 
and the magnetization does not saturate at the highest available field, i.e., 5.5 T. ii) there is 
an important increase of the magnetization at low temperature, iii) the thermal behavior of 
the magnetization at 5.5 T is such that the smaller is the mean diameter of the particle the 
faster is the increase of the magnetization at very low temperature. 
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From the theoretical point of view, the situation involving both surface effects and dipolar 
interactions has never been considered so far mainly because of its tremendous complexities, 
and also because one has first to understand these two effects separately. Needless to say 
that, already at the static level, no exact analytical treatment of any kind is ever possible 
even in the one-spin approximation, i.e., ignoring the internal structure of the particles and 
thereby surface effects. Only numerical approaches such as the Monte Carlo technique can 
relieve some of this frustration. Indeed, applications of this technique to the case of Ising 
dipoles can be found in Ref The same technique has been used in Ref. |3] to study 
hysteretic properties of monodisperse assemblies of nanoparticles with the more realistic 
Heisenberg spin model, in the one-spin approximation where each particle carries a net 
magnetic moment. In Ref. the Landau-Lifshitz thermodynamic perturbation theory 
is used to tackle the case of weakly dipolar-interacting monodisperse assemblies of magnetic 
moments with uniformly or randomly distributed anisotropy axes. The authors studied the 
influence of DDI on the susceptibility and specific heat of the assembly. 

In the present work, we use the same approach as in Ref. Q with the objective to study 
the effect of anisotropy and (weak) dipolar interactions on the field and temperature behav- 
ior of the magnetization of a monodisperse and polydisperse assembly of magnetic moments. 
For this purpose, we consider an assembly of magnetic moments whose magnitudes are dis- 
tributed according to a Gaussian or (the more often observed) lognormal law. The anisotropy 
is taken as uniaxial and either textured along some reference axis or randomly distributed. 
The statistical average of the assembly magnetization is obtained, for weak DDI, using the 
thermodynamic perturbation theory, as in Ref. ^1, but here the magnetic field is explicitly 
included in the assembly magnetic energy. The low field regime, dealt with in Ref. is 
generalized here so as to take account of polydispersity and DDI. In high fields, the mag- 
netization as a function of temperature and field is computed using the steepest-descent 
approximation. In the general range of temperature, field, and anisotropy, the magneti- 
zation of a non-interacting assembly is computed exactly by numerical integration of the 
single-moment (free) partition function or using the Monte-Carlo simulation technique. For 
interacting assemblies, we use the Monte Carlo technique. 

One of the objectives is to provide ready-to- use (semi) analytical formulae for the field and 
temperature dependence of the assembly magnetization that take into account moment and 
easy axes distributions, and weak dipolar interactions. Moreover, we investigate the effect 
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of anisotropy and DDI and discuss the validity of the Langevin law, for both textured and 
random-anisotropy, which is invariably used in the literature to interpret the magnetization 
measurements on nanoparticle assemblies. The present work is also an extension of the study 
in Ref. j^] where anisotropy was ignored. 

The layout of this paper is as follows: Section |H] defines the energy and notation. The 
first new results of the present work appear in section ITTT1 where we use perturbation theory 
to derive an analytical expression for the magnetization taking account of DDI. Then, we 
give approximate expressions in the limiting cases of low and high field regimes, for both a 
free and interacting particle, and for both monodisperse and polydisperse assemblies. We 
compare these expressions with the exact numerical calculation, and also discuss the effect 
of anisotropy. This section ends with a Monte Carlo calculation and discussion of the cor- 
responding results of more realistic assemblies of magnetic nanoparticles, namely assemblies 
with the lognormal distribution for the magnetic moments derived from experimental data 
of nanoparticle assemblies with mean diameter D m = 3 and 7 nm. The last section is a 
conclusion and statement of a few open problems to be dealt with in future investigations. 



Consider an assembly of magnetic moments rrij = m^Sj, % = 1, . . . , Af of magnitude rrii 
and direction Sj, with |sj| = 1. The magnitude of the magnetic moment nij is then defined in 
terms of the Bohr magneton i.e., rrii = riifiB, and the numbers rij are either all equal for 
monodisperse assemblies or chosen according to some distribution, the so-called polydisperse 
assemblies. Each magnetic moment will have a uniaxial easy axis ej, and for an assembly 
these may be either all directed along some reference axis leading to a textured assembly, or 
randomly distributed. The latter case will be referred to as random anisotropy. Hence, the 
energy of a magnetic moment rrij interacting with all the others via DDI, in the magnetic 
field H = He h , reads [after multiplying by — (3 = —l/k B T], 
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is the vector joining the sites i and j and whose magnitude is measured in units of a, 
a characteristic length on the lattice to be evaluated later on. Since we are considering 
assemblies with moment instead of volume distribution, the volume V; in Q is rewritten 
in terms of rii via the saturation magnetization of the material per unit volume M s , i.e., 
Vi = rrii/M s = (nB/M s )rii. For convenience we introduce the dimensionless parameters 

_ n B H fi B K \Xq\i\ 

x ~ k B T' ° M s h B r u wa^t' [ ) 

and define x% = xrii,<Ji = arii. Note that <jj = KVi/(k B T) is the commonly used notation 
for the reduced anisotropy-barrier height of the particle i. Therefore, we rewrite (JIJ as 

Si = Vi* ■ e h + f • e,) 2 + & £ 3(s, ■ e tJ )( Sj ■ e tJ ) - ■ s, (4) 



= ,f 3(sj -e^Csj -ejj) 

— °i ^Z,d2_^ n i' L j 3 



Sj * Sj 



r- ■ 

j<i % 3 



where E^ is the free particle energy. In what follows we also occasionally use the dimen- 
sionless magnetic moment vector Sj = rijSj, and the DDI term is rewritten as a quadratic 
form in Sj 

^MSr^nSyc,,) S,-S, ^ ^ S; £ (5) 
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j<i % 3 j<i j<i 



where we have introduced the DDI tensor 



Vij = {3e i:j e i:j - 1) . (6) 



r ■ 
i.i 



III. MAGNETIZATION 



A. General formulation 



In order to calculate the thermal-equilibrium average of any observable 0{si, . . . , Sj\f) we 
have to average over each particle's moment direction s$, the direction of its easy axis e^, 
and the magnitude m, of its moment (or equivalently n,). The average of O with respect to 
spatial orientations of all spins is 

(0) = ^jvQe £ 0. (7) 

where 

DQe £ , 
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with T>Vt = Y\i d^i — Ylid 2 Si/2TT, and £ = Yli^i- I n particular, the magnetization compo- 
nent along the field taken, for instance, along e z , of a single particle i is given by 

K> = i / Vne £ mt = m,( a j f ), (8) 

which is a function of the easy-axis direction e iy rij and the parameters a, x, ^ (or equivalently 
K, H, T, Next, we infer the magnetization of the assembly per particle as 

1 f d 2 e ^ 

where w(rii) is some distribution of the Bohr magneton numbers ri{. 

Therefore one first has to compute the magnetization of a single particle defined in (JHJ). 
Analytically, this can only be done in the case of weak DDI using thermodynamic pertur- 
bation theory 3|. This operates by expanding the Boltzmann distribution V = Z^ 1 exp(£) 
in powers of the interaction parameter that is 

V = V (l + £ d F 1 + \e d F2 + ...), 

where 

Vo = ^ m ={[vi ) (10) 
Z ° i=i 

is the Boltzmann distribution of the non-interacting (free) ensemble, and 

■W / /■ J2 \ ^ 

*=n(/^)=n* 

with Zq being the directional partition function of the i th free particle. For the system 
considered here, F\ and F 2 are some quadratic, respectively quart ic, functional of Sj. 

Therefore, the calculation of the average of an observable O is reduced to the calculation of 
averages with respect to the distribution Vo of low powers of the spin variables. Henceforth, 
the average with respect to Vo will be denoted by (.) . 

Consequently, to second order in the interaction parameter the average of any physical 
observable O reads, 

(O) ~ (O) + e,A (1) + ^JA {2) + 0(e) (11) 

with, 

' AW = (OGi) - (O) (G l ) , 

(12) 

A^ = (OG 2 ) -(O) (G 2 ) -2(G 1 ) AW, 
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where, 

G 2 = XI $ J + Yl ®ij®kiqik:jiqil:jk, 

i>j i>j k>l 

with q ik .ji annihilating terms containing duplicate pairs: q^ji = |(2 — 5ik — 5ji)(l+5ik)(l+5ji) 
17 1 . Note that in contrast with the situation in Ref . \v\. the l st -order averages (O) do not 
vanish here due to the presence of the external field. 

The magnetization of a single particle interacting with all other particles in the assembly 
can be written to second order in ^ as in Eq. (Jll|) . for the observable Sf, that is the 
(dimensionless) magnetization in the direction of the field taken along the z axis, 

(^)^(^ 2 )o + ^A (1) + ^A (2) , (13) 

where A^ 1 ) and A^ 2 ) are obtained from Eq. (fE?|) by setting O to Sf. 

Now, the calculation of A^ and A^ 2 ^ involves that of averages of products of Sf whose 
order ranges from 1 to 5. Introducing the notation 



(sf) , K 



d(Sf) 
dx 



the average of an arbitrary degree of Sf is expressed by the function hi and its n -order 



derivatives bf^ [2l| . Restricting ourselves to l si -order in and thereby to 3 rd -order averages, 
we have 

(SfSfa = bibj + bfa, (14) 
(SfS*S z k ) Q = bibjbk + b'tb^k + bfaSij + hb'^jk + bffSijSjk. 

Next, noting that i) since the field is applied along the z axis the average of the x, y compo- 
nents vanishes, ii) DDI only involve pairs of distinct indices, say i, j so that Sy = 0, Eq. (fT^j) 
leads to the following expression for the magnetization of an interacting assembly (to first 
order in 

(sn * (sf)o+t d jt,( s k)° A K d< !! f>0 > ( is ) 



dx 

k=l 



with 



[3(e h ■ e kl f - 1] 
s±ki — ^3 — • L>ki ■ eh- 
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As was discussed in Ref.Qand confirmed in section lTlI Cl below. for non spherical systems, 
the corrections to the magnetization are largely dominated by the first order contribution to 
the DDI. 

Upon examining Eq. (fTK|) one sees that the magnetization of an interacting particle is 
written in terms of the magnetization of the free particle and its derivatives. In order to 
render this expression more explicit and thereby more useful from a practical point of view, 
we have to consider some limiting cases where analytical expressions can be derived. So, 
next we compute the free-particle magnetization {S%)o and its derivative with respect to x. 



B. Limiting cases for the free-particle magnetization: effect of anisotropy 

The free-particle (reduced) magnetization is given by 

, „ 1 f tt / d 2 Sk - 1 f d 2 Si F (o) 

This can be computed exactly by numerical integration upon changing to spherical coordi- 
nates. More precisely, using without the DDI term we rewrite (fTo^) as follows 

(s*) = — / dtp due^ y2{u)+x > u Ui (17) 
zq Jo J -i 

where 

y(u) = upi + Vl - u 2 ^Jl - p 2 cos(p, z (mi,ip,a,x) = J dtp J due~i" y2 ^ +XiU , 

with pi = eh ■ Bi. Eq. (fTTj) is valid for all values of a and x (or K, H, T, and Vi). However, 
numerical integration is very time consuming, and knowing that we have to do this for all 
particles and then average over the two distributions, renders this expression of little practical 
interest. Instead, one may derive sensible analytical expressions in the relevant limiting cases 
such as low and high field regimes. The low field case is dealt with perturbatively while the 
high field case is treated using the steepest-descent approximation. 

Let us now give approximate analytical expressions for the magnetization of a particle in 
the assembly in these two limiting cases. 
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1. Low field 

The low-field expansion obtained in reads (upon introducing the particle index i) 



lf 1 + 2S i2 P i2 7 + 70(S i2 P i2 ) 2 + 40S i2 P l2 - 12S i4 Pi 



I/O 



where, 



SilicTi 



315 



(i-W my/a , ^ x 



*4 3 



(19) 



C7i > 1, 



and P*; = Pu{Bh ■ e^) are the Legendre polynomials. For a textured assembly, all of the 
angular functions Pu turn into unity, while for randomly distributed easy axes e, we have 
Pi2 = Pi4 = 0, Pf 2 = 1/5. Hence, the magnetization in low field becomes 



{ l + 23 a 7 + 7051 + 405*2 -125a 3 . , , 
Xi — — X; , textured 



315 



(20) 



I 3 



45 Xi ' 



random. 



It is obvious from these expressions that the magnetization is larger if the anisotropy is 
textured along the applied field. 



2. High field: steepest- descent approximation 



The single-particle partition function Z l Q in Eq. ()16j) reads 



Zl oc / dSi5(ai - l)e 4 



dsj5(s^ — l)e 2 



(21) 



Temporarily dropping the particle index i for simplicity, Z % Q can be rewritten using the 
Hubbard-Stratonovich transformation [22] which consists in introducing an auxiliary (vector) 
field £ and using the Gaussian integration formula 



exp 



° ( \ 
— (s.e) 

2 V ; 




— / duexp( u 2 + a(s.e)u), 

2n 2 



(22) 



so that Eq. (j21|) can be rewritten as 



Z = Const. J due - ** 3 J dsS(s 2 - l)e s ^ = Const, j due~ ps{u ^\ 



(23) 
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where £ is the (auxiliary) effective field acting on each magnetic moment, 

£ = x + true, 



(24) 



and is the (vector) sum of the applied magnetic field and the anisotropy field; S(u,£) is the 
effective action given by 



The variable u appearing in Eq. 
which thereby reads 



ds5(s 2 



lie 



Const. + ln 



sinh(£) 



(25) 



is determined through the minimization of S(u, £), and 



u 



Co 



-Co • e, 



(26) 



where -B(x) = coth(x) — •= is the Langevin function. This is a transcendental equation for 
u Q leading to a transcendental equation for £ , 

5(6)) 



£o = x + a^^(e -e)e. (27) 
so 

23j | to compute the par- 



Next, we use the steepest-descent (or saddle point) approximation 
tition function and the action of the particle in a high magnetic field, i.e., x ^> 1. This 
consists in Taylor expanding the action f!25D around uq, which then may be rewritten as 



Co 



1 - 



(Co • e) 



(28) 



which can be further expanded. Since, in terms of (« ,(o), the free energy and the action 
are equal, the particle magnetic moments is given by 

~d(f3S)~ 



<9x 



(29) 



(30) 



and along the field direction, i.e., e z = e h = x/x, we get 

(s z ) = (s).e h = ^-(^.e h ). 
so 

In order to compute (s z ), we must solve (}2Tj) for £ , and this can be done perturbatively 
assuming that the applied field is large and proceeding by successive expansions in 1/x and 
replacing several times uq and Co by their expressions (J26|) and P7jl . respectively. Note that 
Co depends on -B(Co) but since the Langevin function rapidly saturates to 1, in the high-field 
case we may take -B(£ ) = B(x) or even £?(Co) = 1 — 1/x. Consequently, we obtain 



£ = x + op 2 



1 + 



-i + Ki-pV 
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Upon using these expressions and reinstating the moment index i, we obtain the approximate 
expression for the free particle's magnetization as an expansion in l/x, 



(sX - 1 - J. + (31) 



where, 



Si = ^ [3p 4 2 - 1 + apiipl - 1)] , Pi = e h - e, t . 



2 

Eq. ({HI} is only valid at high fields but at all temperatures, though it yields a better 
approximation at low temperatures. More precisely, it is valid for fields H larger than 

H = ^x ■ x =i±^5 A = 1-4S. 

/iBUi 2 

Expression (|ST|) is much easier to handle than (fTTj). and it yields the same result in its range 
of validity, i.e., for x > x m \ n . 

For a textured and random-anisotropy assembly, we infer from Eq. (J31|) 

f 1 1 °i 
1 1 — 25 textured 



(32) 



random. 

a;,- 15 27 



The validity of the asymptotic low field and high field expressions of the magnetization 
given in (J2U|) and respectively, is checked by comparing the latter to the exact nu- 

merical results obtained from Eq. (|17j) . For example, Fig. ^ shows such a comparison for a 
monodisperse assembly and random anisotropy. It is clear that the asymptotic expressions 
are good enough. Moreover, even at a relatively high temperature (here small a), for which 
the steep est- descent approximation is expected to work worst because x becomes small, 
expression renders a good approximation. 



3. Comparison with the exact calculation and Langevin's function 

We can immediately infer some results from the low and high field expressions derived 
above. Indeed, considering that the Langevin function £(xi) = cothxj — 1/xj expands to 
~ Xi/3 — at low fields and to ~ 1 — 1/xj at high fields, Eqs. (J2T1|) . (j22I) readily imply 

that: 
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Reduced magnetization (per particle) 




FIG. 1: Exact (numerical) calculation versus approximate expressions of the reduced magnetization 
of a monodisperse assembly with random anisotropy. 

• both low-field and high-field curves fall onto the Langevin curve if the assembly is 
isotropic, i.e., crj = 0. 

• the magnetization of a textured assembly falls above the Langevin curve, while that 
of an assembly with random anisotropy falls below Langevin's curve, and the larger is 
the anisotropy (hence <7j), the larger is the deviation. 

To confirm these results, we numerically compute the magnetization (fTTj) for monodisperse 
and polydisperse assemblies of magnetic moments, with random anisotropy. In the second 
case we considered a simple Gaussian (or normal) distribution. The results are presented in 
Fig. 121 This clearly shows that as a increases, or equivalently at a fixed temperature and 
increasing anisotropy constant, the magnetization drops. Indeed, stronger anisotropy implies 
that it is more favorable for the magnetic moments to align along their randomly oriented 
easy axes, and so the Zeeman energy is not sufficient to align them along the field direction, 
as is clearly shown by Eq. (jlOJ), (jHU)- This holds for both monodisperse and polydisperse 
assemblies. Note also that the zero or very weak anisotropy curves coincide with the Langevin 
function £(£) = coth(£) — l/£, which simply confirms the fact that Langevin's law is only 
rigorously valid in the absence of anisotropy or at high temperature, or more precisely in 
the superparamagnetic regime. On the contrary, for a textured assembly [whose results are 
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FIG. 2: Reduced magnetization (per particle) of an assembly with randomly distributed anisotropy 
axes. Left: for a monodisperse assembly. Right: for a polydisperse assembly with Gaussian 
distribution for the magnetic moments. 



not shown here] all easy axes are parallel and obviously stronger anisotropy leads to larger 
magnetization. 

We have also investigated the effect of anisotropy on the magnetization of some typical 
nanoparticle assemblies with the most often observed volume (or n) distribution in samples 
of magnetic nanoparticles, namely the lognormal distribution with parameters fi, 5, 

2" 



w(n) 



1 



cxp 



1 /In n — /j, 

2 V 5 



(33) 



We consider assemblies of, e.g., cobalt or maghemite, nanoparticles with mean diameter 
D m = 3 and 7 nm. We simulate the assembly as a collection of magnetic moments randomly 
assigned to the sites of a regular simple cubic lattice. The moment magnitude distribution 



is taken from experiments 



24j upon converting the volume or diameter to the corresponding 



number of Bohr magnetons n. The magnetization of these assemblies is computed using the 



standard equilibrium Monte Carlo technique 



14. 



at arbitrary temperature, applied field, 



and DDI parameter For a non-interacting assembly the field behavior of the magneti- 
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Reduced magnetization (per particle) 




T = 5K 
Polydisperse assembly 
Random anisotropy 



40 X 



50 



FIG. 3: a) Reduced magnetization (per particle) of an assembly of M = 1024 with lognormal- 
distributed magnetic moments with mean diameter D m = 7 nm and randomly-distributed easy 
axes as obtained from Monte Carlo calculations for different anisotropy values. x m = n m HBH JksT ', 
where n m is the mean number of Bohr magnetons for this assembly, b) Langevin function together 
with the Monte Carlo results for D m = 3 and 7 nm. 

zation, for different values of the parameter a is shown in Figs. EJ In Fig. we observe 
that, in the high field regime, the conclusions drawn from Fig.|2]are confirmed in the present 
case too. That is, the higher is a, the lower is the magnetization. On the other hand, 
in low fields this is not globally so because the competition between Zeeman, thermal and 
anisotropy contributions to the energy, for this distribution, results in a crossing between 
the various magnetization curves, as has been observed, e.g., for maghemite particles [see 
(Fig. 10) and 12 1]. In fact, this situation is reminiscent of the two phases [blocked and 
superparamagnetic] exhibited by the zero-field-cooled magnetization, and separated by the 
temperature T max at the peak. More precisely, at a given applied field H, if x is sufficiently 
large, T is smaller than T max , and thus the larger is the anisotropy and/or particle volume, 
the higher is the energy barrier, and the smaller is the magnetization, because most of the 
particles remain in their blocked states with almost randomly oriented moments. In the 
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opposite situation, small x corresponds to high temperature, and in this case anisotropy 
has a negligible effect, so that a higher volume corresponds to a higher magnetic moment 
(or n) and thereby to a higher Zeeman energy, which is necessary to take over the strong 
thermal fluctuations. Higher Zeeman energy, of course, builds larger magnetization. Equally 
important is the observation of the related effect, also observed experimentally in maghemite 
particles [see e.g., Fig. 10 of Ref. 0], that the magnetization as a function of field [see Fig. EK| 
has a much larger slope for small x than for large x. Furthermore, we see that there is a large 
deviation from the Langevin law, due to several parameters, ignored by the latter, especially 
anisotropy. In addition, we see in Fig. Eb that, as was observed earlier, the larger is the 
mean diameter of the assembly, the larger is a, and thereby the larger is the deviation from 
the Langevin curve. 



C. Effect of dipolar interactions 

Now, we derive the expressions analogous to Eqs. (f2UJ) . (f3*2j) for a weakly interacting 
polydisperse assembly. We only do this in the case of random anisotropy. Accordingly, in- 
serting the low and high field expansions (J2U|) . in Eq. (fTB^) leads to analytical expressions 
for the magnetization of a weakly interacting assembly as a function of field, temperature, 
anisotropy, and the DDI parameter The assembly (reduced) magnetization per particle 
is defined as (s z ) ass = l/A/"X}i=i ( s z)i- In the case of randomly distributed easy axes we 
obtain, 



1 + l^ 1 - 2 ) 



3 



4~ 

A 3 + -i d A b 



{xf 
45 



(x) 



15 



(xy 



2(^) 2 c (i,ii 
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£(0,0) 



(X) 



3 ' 



low field 



high field, 



(34) 

where we have defined [see Eq. for notation] ^ = ^ (n) 2 , (x) = (n) x, (cr) = (n) a, with 
(n) = l/J\fYli=i n ii an d the (scaled) constants 



N . rr*., . (n) a+b ' 3 *f^{n) 



i M 3 



3 ! 



A 5 - 
5 Af 



1 ^ ^nfAjrf 
Tf 2. (n> 5 • ^ 



on = 1 + 2Sf 2 with S i2 defined in dHJ). 
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Note that there are three types of contributions: There are pure anisotropy terms, pure 
DDI terms, and mixed terms. It is readily seen that in the absence of anisotropy and DDI, 
i.e., for crj = 0,^ = 0, and from Eq. (fTH|) ai = 1, expressions (pMj) simplify back, in the case 
of a monodisperse assembly with rii = n, to the expansions of the Langevin function in low 
and high field regime given at the beginning of section IIII B 3| or as can be inferred from 
Eqs. (JUJ), (jH2I)- Let us now discuss the constants C (ct ' 6) , A 3 , A 5 . 
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FIG. 4: Some of the constants of Eq. versus the vertical size Z of the cubic lattice 

10 x 10 x Z. For Z < 10, the lattice is oblate, for Z > 10 it is prolate, and for Z = 10 it is 
cubic. These results are for a cobalt assembly with D m = 3 nm (thick lines), and 7 nm (symbols), 
corresponding to n m = (M s V m / ' Hb) — 2172 and n m ~ 27595, respectively, where we have used 
V m = ttD^/6, and M s ~ 1425 x 10 3 J/T/m 3 . 

First, in the continuum limit the lattice sum C^ ' ^ = z~2ij=i i^j A%j becomes C^ ' ^ = 

47r(l/3 — A 2 ), for a simple cubic lattice, with A 2 being the demagnetizing factor along the z 

axis. For instance, for a box with semi-axes a = b = 5, c = 10, we have 

abc f°° ds 



X: 



0.174, 



2 Jo (c 2 + s)v/(a 2 + s)(6 2 + s)(c 2 + s) 
so that = 47r(l/3 — X z ) ~ 2. Note that in the monodisperse case we haveC^ oc . 

In the polydisperse case, using the assemblies of Fig. in Fig. |U we plot these constants 
as functions of the size Z along the field direction for a box-shaped lattice with size M = 



17 



X x Y x Z, and X = K = 10, for the mean diameter of 3 and 7 nm. It is clear from 
Fig. E| that these constants, and thereby the corresponding DDI terms in Eq. ([34)1. are shape 
dependent. This is no surprise knowing that the long range DDI lead to shape dependence 
of the physical quantities, and in particular the magnetization. On the other hand, we see 
that these constants are negative for the oblate system, and positive for the prolate, which 
implies that the DDI suppress the assembly magnetization in the former case and enhance it 
in the latter. Moreover, it is also seen that for cubic systems all the constants C( a > fc ) vanish, 
which means that the DDI do not contribute to the magnetization in this case, and thus 
the deviations from the Langevin behavior are caused only by anisotropy. Note that all the 
scaled constants C^' 6 -* are almost independent of the assembly mean diameter. In particular, 
this is trivial for C^ ' -*. 

In order to study the constants A 3 , A 5 in a similar way, we have to make assumptions 
about the intensity of the anisotropy, since these constants contain the parameter «j. From 
Eq. ()19)) we infer that in the absence of anisotropy, i.e., Oi = 0, «j = 1, while for strong 
anisotropy we may approximate Sa to 1, and hence to 3, so that A 3 oc (l/AA) X}i n ?/ ( n ) 3 
and A§ oc C*- 2,3 ^ in both limits of anisotropy. For D m = 3, 7 nm, A3 ~ 2, 6. In the continuum 
limit A 3 tends to exp(35 2 ), where 5 is the standard deviation of the distribution ()33|) . A 5 
shows the same behavior as C^ 1 ' 2 ' but with bigger change with D m . It is well known from 
other areas of physics that the calculation of such high-order moments (or "cumulants") 
requires more precision because they present more statistical fluctuations with the lattice 
size. 

In Fig. 03 we plot the Langevin function (full line) and the Monte Carlo results (symbols) 
for the magnetization of an interacting assembly of (TV" = 10 x 10 x 5) lognormal-distributed 
moments, with random anisotropy, and for different values of the inter-particle distance. 
Here we use the same assemblies as in Fig. El The intensity of DDI, or equivalently the value 
of is varied by varying the lattice parameter a entering ^ [see Eq. (JHJ)]. More precisely, 
the parameter a is taken as a real number times the mean diameter D m of the assembly, i.e., 
a = k x D m . Thus, large values of k correspond to an isotropically inflated lattice with large 
distances between the magnetic moments, and thereby weak DDI. 

These results, obtained for an assembly on a simple cubic lattice, do confirm that DDI 
suppress the magnetization. Indeed, we recall that it was shown by Luttinger and Tisza 
|28| [see also the more recent work jlfjj using the Monte Carlo technique] that the ground 



18 



state of a simple cubic lattice of dipoles is antiferro magnetic, while that of a face-centered 
cubic lattice is ferromagnetic. It is also seen that these curves deviate from the Langevin 
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FIG. 5: Reduced magnetization (per particle) of an interacting assembly of N = 10 x 10 x 5 
lognormal-distributed magnetic moments with mean diameter Dm = 7nm and random anisotropy. 
Monte-Carlo in symbols and analytical expressions (|34j) in lines. The parameters £ and k are 
defined in the text. 



law. However, we emphasize that the deviations induced by DDI are much smaller than 
those induced by anisotropy, as already discussed earlier [see also Ref. Q for a related 
discussion of the effect of the system shape on the magnetic susceptibility of a monodisperse 
assembly]. In the inset of Fig. El the same results are magnified by plotting them in function 
of £ = {(%) /id) x oc fiBH/(fi 2 B /a 3 ), i.e., the ratio of Zeeman energy to the DDI energy, 
which also makes it possible to distinctively plot the analytical expressions (jSljl for low field. 
In the case of high field only one curve (k = 2, i.e., relatively strong DDI) is presented since 
for the other values of DDI parameter k, the steepest-descent approximation is valid for 
much higher values of (. Note also that in function of the parameter ( the tendency with 
increasing DDI strength is reversed. 
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IV. CONCLUDING REMARKS AND OPEN PROBLEMS 

We have provided simple approximate analytical expressions for the magnetization of a 
weakly interacting polydisperse assembly of magnetic moments with randomly distributed 
easy axes, in both low and high field regimes. These expressions have been checked against 
extact calculations using either numerical integration of the partition function or Monte 
Carlo simulations. We have also computed the magnetization of such systems for an arbi- 
trary inter-particle separation, or equivalently arbitrary intensity of dipole-dipole interac- 
tions, using the Monte Carlo technique. However, in this case only assembly of limited sizes 
could be dealt with, as the calculation speed is reduced to a crawl by the long-ranged DDI. 
We have investigated the deviations caused by random anisotropy, DDI, and polydispersity, 
from the Langevin law that is commonly invariably used in the literature to describe the 
magnetization of real materials. We have also found that more realistic assemblies with a 
lognormal volume distribution, render a magnetization that exhibits two different regimes 
as a function of the applied field, with different variation slopes, as has been observed in 
experiments on maghemite particles. Moreover, as a byproduct, we find that the magneti- 
zation of an assembly of nanoparticles in the one-spin approximation considered here, i.e., 
with each particle represented by a macro-spin, does saturate in high fields. This suggests 
that the magnetization non-saturation observed in experiments on small particles is most 
likely due to some intrinsic properties of the particles, such as surface effects, as has been 
argued in many publications [see e.g., 0, 0| and many references therein]. 

For future investigation, we intend to apply the Fast-Fourier- Transform technique to speed 
up the Monte Carlo calculations especially for interacting assemblies of more reasonable sizes, 
and take account, inter alia, of random spatial distributions of the particles on the lattice. 
Using the kinetic Monte Carlo technique, we also intend to investigate the disappearance 
of the maximum in T max (H) as the concentration of the assembly is increased, and also the 
appearance of the spin-glass like state at low temperature. For the latter purpose, we will 
most likely have to tackle the problem of an assembly of multi-spin particles. 
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